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Abstract 

The explicit dispersion equation for a one-dimensional periodic structure with alternative layers 
of left-handed material (LHM) and right-handed material (RHM) is given and analyzed. Some 
abnormal phenomena such as spurious modes with complex frequencies, discrete modes and photon 
tunnelling modes are observed in the band structure. The existence of spurious modes with complex 
frequencies is a common problem in the calculation of the band structure for such a photonic crystal. 
Physical explanation and significance are given for the discrete modes (with real values of wave 
number) and photon tunnelling propagation modes (with imaginary wave numbers in a limited 
region) . 

PACS numbers: 78.20.Ci, 42.70.Qs, 78.20.-e, 03.65. Ge, 
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Left-handed materials (LHM) with negative permittivity and negative permeabihty, 
which were first suggested in have attracted much attention recently 0, ^ |[ thanks 
to the experimental realization of such materials and the debate on the use of a LHM 
slab as a perfect lens to focus both propagating waves and evanescent waves 0, • 

Multi-layered structures containing LHM have been investigated through calculating the 



transmittance or refiectance of the structure |TU|. Effects of photon tunnelling and Bragg 
diffraction have been observed in these works. In the present paper, we consider a one- 
dimensional (ID) photonic crystal formed by alternative layers of LHM and RHM. The 
band structure is investigated by analyzing the explicit dispersion equation. Some abnormal 
phenomena such as spurious modes with complex frequencies, discrete modes and photon 
tunnelling modes are found and explained in the present paper. 

Consider a ID periodic structure with alternating layers (ei,/ii) and (e2,/i2) as shown in 
Fig. ||. di and d2 are the widths of the two inclusion layers respectively, and a = di + d2 
is the period. Note that the conservation of energy requires that ej/ij > 0, i = 1, 2. The 
corresponding refractive index is given by rii = ±y/eifii (negative sign for LHM). We consider 
an off-line propagation of monochromatic electromagnetic field (with time-dependence e"""^*) 
in a periodic structure with off-line wave vector j3 along the x axis. For the E-polarization 
case, one has 



Ely = e'P^[eik^^ Ae-'^ 



Hi^ = ^e'^''{e'''^^ - Ae-^^^^) in region 1, (1) 



and 



= e'^'^iBe'^^^ + Ce''^^^) 

= ds^f^iUxi^Be'^^"^ - Ce^*^22) in region 2, (2) 
= ^e^^^(5e*^22 _^ Ce-'''^') 



where ki is the component of the wave vector along the axis in region i {i = 1,2), i.e., 
kf = u'^eiHi/c^ — P"^. Here c is the speed of light in vacuum. 

The tangential electric and magnetic fields should be continuous at z = 0, i.e.. 



(3) 
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To obtain the dispersion relation for this ID photonic crystal, we need to use the following 
periodic conditions according to the Bloch theorem, 

-- d2) = Eiy{z = -di)e"^- 

(4) 

= d2) = Hix{z = -di)e'^% 

where q is in the first Brillouin zone — ^ < g < ^ . Substituting systems (1) and (2) into 
systems (3) and (4), we obtain the following dispersion relation for the E-polarization case, 

cos(A;ic?i) cos{nkid2) — : sin(fci(ii) sm.{nkid2) = cos(ga), (5) 




where h = ^2/^1 = \Jn'^{k'l + — /i = /i2//^i, and n = n2/ni is the index ratio of 

the two media. 

Similarly, we can derive the following dispersion relation for the H-polarization case. 
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cos{kidi) cos{nkid2) — — sm{kidi) sin{nkid2) = cos(ga), (6) 

where e = e2/ei. 

From analytical equations (5) and (6), we can find out how ki depends on q. The 
corresponding dispersion relation between uj and q can then be obtained according to u;^ = 
c^{kl + /5^)/(ei/ii) for each fixed /? (/5 = is associated to the in-line propagation, i.e, the 
case when the propagation direction is along the normal of the material interfaces). 

When the index ratio n < 0, the solution to the dispersion equation and the associated 
band structure have quite different behaviors as compared to those for the usual case when 
n > 0. For simplicity, we first illustrate this for the case when /3 = 0. In this case. 



2hfj. ^ dispersion equations for the E-polarization and the H-polarization 

have the same form (as expected). 

Spurious modes with complex u. When one of the constitutive media is of LHM 
(i.e., n < 0), the analytical dispersion equation may have some complex solutions for ki. 
These complex values of ki lead to complex u since u"^ = c^{k\ + [S"^) / {tijii) . Fig. || gives 
the in-line (i.e., /5 = 0) band structure when n = — 2.5,yU = —2,di = ^2 = 0.5a. The real 
part of uj varies little as q varies from — vr/a to vr/a. The problem of complex u will also 
appear when /3 7^ (i.e., the case of off-line propagation). Obviously, these complex u have 
no physical significance since physically the frequency u must have a real value. 

One would meet a similar problem of complex u in the calculation of the band structure 
using other methods such as the plane wave expansion method (one needs to remove these 



spurious complex uo in the band structure). In the plane wave expansion method, the uj 
values are determined from the following eigenvalue equation, x (^V x H) = ^H. The 
existence of complex uj indicates that the operator x (^V) is not a Hermit operator or 
has no equivalent Hermit operator. One can also expect that these complex uj (associated 
with an exponential increasing factor of time) will make the finite-difference time domain 
(FDTD) method (with a periodic boundary condition) divergent in the calculation of the 
band structure of a photonic crystal consisting of alternate RHM and LHM inclusions (with 



constant permittivity and permeability) |jTT], 12, 13 1. We have verified numerically with the 



plane wave expansion method and the FDTD method the existence of the spurious modes 
with complex uj. 

The fields associated with the spurious complex oj satisfy Maxwell's equations, the bound- 
ary condition and the enforced periodic condition. Therefore, an additional requirement that 
the frequency uo must be real should be applied to make the solution significant in physics. 

In the rest of the paper, we consider only situations when uo is real. 

When /3 = 0, a complex solution of ki always leads to a spurious complex uj. However, 
when /3 7^ some complex solutions of ki may correspond to real uj and thus have physical 
significance. Since uj"^ = c^{kf + j3'^)/{eifii), a complex solution of k^ = ki^ + jku (with 
kii 7^ 0) should satisfies the following conditions in order to make uj real, 

A;,^ = 0, kl<^\ (7) 

In other words, an imaginary solution of ki can have physical significance when /5 7^ (i.e., 
the off-line propagation case). Imaginary solutions of ki correspond to evanescent waves. 
In a periodic structure consisting of alternate RHM and LHM layers, the special property 
that the evanescent waves (in both directions) decrease in the RHM layer and increase in 
the LHM layer |^ or vice versa permits the evanescent waves to propagate in the periodic 
structure (i.e., a propagation mode, for which the energy density does not decay over a 
period; we will call it a photon tunnelling mode and will discuss it in more details below). 
In the periodic structure composed solely by RHM or LHM, complex solutions of ki will not 
occur since evanescent waves can not form a propagation mode in these structures. 

As discussed before, imaginary solutions of ki may have physical significance when /3 7^ 0. 
This will be illustrated by the abnormal photon tunnelling modes in the band structure for 
a photonic crystal consisting of alternate RHM and LHM layers. Another abnormal (and 



interesting) phenomenon is that some discrete modes (with real values of wave number) are 
also found in the band structure. 

To simplify the analysis, we consider a special case when n = —1 for the E-polarization. 
Then n = —1 (independent of (3) and Eq. (^ is simplified to 

1 + u2 

cos{kidi) cos{kid2) sin{kidi) sm{kid2) = cos(ga). (8) 

^ IJu 

One notices that the above equation for ki is independent of /3.We discuss situations when 
yU 7^ — 1 (so that the RHM and LHM layers have different values of impedance). It is 
particularly interesting when the two inclusion layers have the same width (i.e., di = d2 = 
0.5a) since in this case only discrete modes and photon tunnelling modes exist (conventional 
propagating band disappears in the band structure). 

Discrete modes. Since cos(fci(ii) cos(/i;i(i2) + ^^j^ sin(fci(ii) sin(A;i(i2) > 1 (when ki is 
real), real solutions of ki to Eq. (8) exist only when g = (notice that —n/a < q < it /a), 
which corresponds to discrete solutions ki = 2NTT/a, N = 0, ±1, ±2, ■ ■ -. Since these discrete 
solutions satisfy the transparent condition of kidi = Ntt, k2d2 = —Ntt (one can easily show 
that each layer is transparent to neighboring layers at these discrete frequencies for a fixed (3), 
the field with this wave vector will surely pass through the periodic structure and becomes 
a propagating mode. On the other hand, fields with any other wave vector (i.e., any other 
frequency for the fixed (3) satisfies the Bragg refiective condition kidi + k2d2 = and can 
not propagate in this periodic structure. Therefore, only those discrete propagation modes 
exist. As an example, the discrete real solutions of ki to Eq. (8) are shown by the circles 
in Fig. ^ (with fi = —2). This abnormal phenomenon (existing for any (3) can be utilized to 
make a very narrow filter (with no sidelope, which is different from any conventional type 
of filters). 

Photon tunnelling modes. From Eq. (8) one sees that only complex ki can be the 
solutions if g 7^ 0. Then according to condition (7) the solutions ki with physical significance 
must be imaginary. The expression on the left-hand side of Eq. (P) always decreases when 
kii increases or decreases away from (a symmetric function of ku with a maximum at 
kii = 0). Therefore, only two conjugate imaginary solutions exist. The dashed lines in Fig. |^ 
show the imaginary solutions of ki to Eq. (8) for an example. When 0) varies in the 
first Brillouin zone, the two conjugate imaginary solutions vary continuously in a region from 
— 0.561jx27r/a to 0.561jx27r/a. Note that this figure is independent of /5 for this special case. 
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For off-line waves with j3'^ > kfj, the imaginary solutions correspond to special propagation 
modes— photon tunneling modes (as discussed before). For this example, evanescent modes 
with kij outside the region [— 0.561j x 27r/a, 0.561j x 27[/a] can not propagate inside the 
periodic structure. From Fig. ^one can also see that ku approaches when q approaches 0. 
Thus, any (3^0 can permit the existence of some photon tunnelling modes since ku can be 
arbitrarily small to satisfy kfj < (3'^. From the above discussion, one sees that for any /? 7^ 
there always exist some evanescent modes which can propagate inside the periodic structure 
while some other evanescent modes can not (one notices that in the case of a perfect lens ^ 
all the evanescent waves can go through the lens). This limitation of the photon tunnelling 
modes is due to the enforced periodic condition. Thus, the problem of square non-integrable 
field mentioned in will not occur in the present situation (this is also true for the next 
example when di 7^ 6/2). 

Both discrete modes and photon tunnelling modes exist when /5 7^ 0. However, when 
/? = (as in the above example), which corresponds to normally incident waves (i.e., the in- 
line propagation case), only discrete modes can exist (this means that only fields at discrete 
frequencies can exist in the structure). 

In the above discussion the case when di = ^2 is considered. Fig. ^(a) shows the solutions 
of ki for an example when di 7^ d2- From this figure one sees that non-discrete real ki 
solutions (corresponding to some continuous bands shown in Fig. ^(b) exist besides some 
discrete modes and photon tunnelling modes. The discrete solutions are at ki = (2.5 + 
5A^) X 2TT/a for q = ±-n/a and at A;i = 5A^ x 2-n/a for g = 0, = 0,±1,±2,. They are 
located in the forbidden band-gaps (see Fig. §(b) for the corresponding band structure for 
a given (3). 

In the above we have discussed only the ideal situations when the periodic structure 
consists of ideal materials with lossless and non-dispersive negative e and negative /i (e and 
are the relative permittivity and permeability of the two inclusion media). In a non-ideal 
situation with material loss or dispersion the abnormal phenomena such as spurious modes 
with complex discrete modes and photon tunnelling modes can also be expected to occur 
for the following reasons. 

If the inclusion media are dispersive, e and /i will be frequency-dependent. Since the 
solutions can be obtained from Eq. (5) in a point-by-point manner, they are correct lo- 
cally (i.e., in a small frequency range where the frequency-dependence of e and /i can be 
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neglected). Thus, the abnormal phenomena can still exist in certain small frequency region. 
Furthermore, the frequency in the band structure is scalable and thus we can scale it to a 
frequency at which the material parameters have desired values in order to observe these 
abnormal phenomena. 

If the LHM inclusion is lossy the relative e or /i has a very small positive imaginary 
part (e.g., 10~^j; as claimed in some papers that LHM must have some loss |^). The 
abnormal phenomena can still occur since the dispersion equation (Eq. (5)) will not give 
any rapid change when a very small imaginary part is added to e or /x. Note that for the 
photon tunnelling modes q should be complex (corresponding to a pseudo-periodic field 
with a damping factor along the propagation direction) when e or /i is complex. Since 
the imaginary parts of these parameters are very small, Eq. (5) is still valid when all the 
parameters (including e, and q) are replaced with their real parts. 

In conclusion, we have found some abnormal phenomena such as spurious modes with 
complex discrete modes and photon tunnelling modes by analyzing the explicit dispersion 
equation for a ID periodic structure with alternating LHM and RHM layers. One will meet 
the problem of complex uo in the calculation of the band structure using other methods such 
as the plane wave expansion method and the FDTD method. The physical significance of 
the discrete modes and photon tunnelling modes has been explained. The discrete modes 
can be utilized to make a very narrow filter with no sidelope. For an off-line wave there 
exist some photon tunnelling propagation modes while some other evanescent modes can 
not propagate inside the periodic structure. 
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FIG. 2: The band structure for the in-line propagation (i.e., /? = 0) with ei = = l,n = 

— 2.5,/i2 = — 2,di = d2 = 0.5a. The circles show the real while the dashed lines indicate the 
real parts of the complex uj. 
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FIG. 3: The solutions (with physical significance) of ki to Eq. (8) with ei = l,/xi = l,n = — 1,/X2 = 
— 2,di = 6,2 = 0.5a. The circles indicate the real solutions, which exist discretely at g = 0. The 
dashed lines show the imaginary solutions. The imaginary solutions are continuously distributed 
in a region from — 0.561j x 27r/a to 0.561j x 27r/a. 



9 




qa/2ji qa/27t 

FIG. 4: Off-line case witli di ^ 3,2- The two inclusion layers have ei = 1, ^ui = 1, n = —l,fi2 = 
—2,di = 0.4a, and d2 = 0.6a. (a) The solutions of ki with physical significance. The circles show 
the real solutions. The real solutions have discrete values (at g = and q = zt-Tr/a) as well as 
continuous bands. The dashed lines indicate the imaginary parts of the imaginary solutions, (b) 
The band structure when /3 = 0.583 x 27r/a. The dashed lines show the frequencies of the photon 
tunnelling modes. The circles indicate the discrete frequencies of the discrete modes. The solid 
lines show the continuous bands. 
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